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Abstract 

Let Aq{a' ,a,k) be the scattering amplitude, corresponding to a local potential 
q{x), X G M^, q{x) = for |x| > a, where a > is a fixed number, a', a £ 5^ are 
unit vectors, S'^ is the unit sphere in M^, a is the direction of the incident wave, 
/c^ > is the energy. We prove that given an arbitrary function /(a') G L^(S'^), 
an arbitrary fixed ceo £ 5'^, an arbitrary fixed k > 0, and an arbitrary small e > 0, 
there exists a potential q{x) G L^{D), where D C ii^ is a bounded domain such 
that 

\\Aq{a ,ao,k) - f{a)\\L2(^s^) < e. (*) 

The potential q, for which (*) holds, is nonunique. We give a method for finding 
q, and a formula for such a q. 

1 Introduction 



If q{x) = for |x| > a, where a > is some number, q{x) G L'^{Ba), Ba = {x : \x\ < a}, 
then the corresponding scattering amplitude Aq{a', a, k) is defined as follows. Let a G 5^ 
be a given unit vector, where 5*^ is the unit sphere. The scattering problem consists in 
finding the scattering solution u{x, a, k), which solves the equation 

[V^ + P - = in M^ (1) 
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and satisfies the following asymptotics: 
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Vector a is called the incident direction, the direction of the incident plane wave e , 
vector a' is the direction in which the incident wave is scattered, and the coefficient 
Aq{a', a, k) is called the scattering amplitude. The properties of the scattering amplitude, 
corresponding to a real-valued, rapidly decaying at infinity q, has been studied in detail 
in the literature. (See, e. g., |21). 

The inverse scattering problem with fixed-energy data consists in finding a potential 
q{x) from the knowledge of Ag{a', a, k) at a fixed A; > and all a', a G S"^. The uniqueness 
of the solution to this probem in the class of real-valued, compactly supported, square- 
integrable potentials was first proved by the author ([4J), who also gave a method for 
recovery of q from the exact fixed-energy data, an error estimate for this method, and 
stability estimates for the recovery problem ([7j) as well as a method for a stable recovery 
of q from noisy fixed-energy data, and an error estimate for this method jS] , jB] , [Z] • Until 
now this is the only known rigorously justified method for recovery of q from noisy fixed- 
energy data. 

In this paper a new problem is studied. Let us assume that the incident direction a 
is fixed, a = ao, and > is also fixed, k = k^ > 0. Denote by Aq{a') := Aq{a', a^, ko) 
the corresponding scattering amplitude, and let /? := a'. Let G L'^{S'^) be an 

arbitrary function. In part of our arguments, in Lemma 3, specifically, it is convenient 
to assume first that the norm of / is sufficiently small norm ||/(/5)||j;,2(52). The role of 
this "smallness" assumption will be clear from our arguments. Under this "smallness" 
assumption we derive an analytical explicit formula for the potential which generates the 
scattering amplitude at a fixed > and a fixed incident direction a with any desired 
accuracy. The "smallness" assumptions will be dropped in Lemma 4. 

Fix an arbitrary small number e > 0. Let D C be a bounded domain. 

The problem (P) is: 

Given e > 0, /(/?) G L'^{S'^), 1 1^2(52), ao G S'^, and ko > 0, does there exist a 
potential q G L'^{D) such that 



and how does one calculate such a q? 

Our basic result consists of an answer to these questions. In jS] the approximation 
problem problem similar to the above was considered for the first time, but the argument 
in [H] does require the "smallness" assumption. 

Theorem 1. For any f G L^(S'^), an arbitrary small e > 0, any fixed a = & S"^, any 
fixed k = ko > 0, and any bounded domain D C M^, there exists a (non-unique) potential 
q{x) G L'^{D), such that © holds. 



||/(/3)-A,(/?)|U2(52)<e, 



(3) 
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To calculate such a potential, we need some auxiliary results. 



Lemma 1. Given an arbitrary f\[3) G L'^{S'^) and an arbitrary small £ > 0, there exists 
an h E L'^{D), such that 



J D 



< e. (4) 



In Lemma 1, and in Lemma 2 below, we do not need the "smallness" assumption. In 
what follows we always assume that a = ao € S*^ and /c = fco > are fixed, and write 
a, in place of and k^. 

The following formula for the scattering amplitude is well known: 

Ag{/3,a,k) = -^ [ e-'^^-''q{x)u{x)dx, (5) 

where u{x) := u{x, a, k) is the scattering solution and the dependence on a and k is not 
shown because a and k are fixed. 
If 

h{x) = q{x)u{x), (6) 

then (jH) is identical to (jHI). 

We now explain how one calculates q from equation © if is given, and how one 
calculates h, satisfying (jH), if / is given. 

Lemma 2. Given f G L^(S'^), one calculates h, satisfying (Q, by the formula: 



i2 I < L 

j-kSi.i^c^y ~ ' (7) 

, i>L 



where g^i^vik) := x^^^ J^(^kx)dx. This integral can be calculated analytically ((T], for- 
mula 8.5.8), and we have assumed that h{x) = for \x\ > 1. 

If h is found such that (jl)) holds, then q can be calculated from equation (jHl), which 
is a nonlinear equation for q, because the scattering solution u depends on q. Under the 
" smallness" assumption we derive an analytical explicit formula for q, namely, formula Q 
below. This formula holds not necessarily under the "smallness" assumption, as follows 
from Lemma 4. 

In Lemma 3 we use the "smallness" assumption for the first time. This assumption 
guarantees that inequality (jH)) holds. This inequality is sufficient for the formula to 
produce an L'^{D) potential q. 

In Lemma 4 we prove that there exists an /i^, a small perturbation of | /i<5| |l2(d) < 
5, which leads by formula ©, with h replaced by /i^, to a potential g^, generating the 
radiation pattern as close to the desired / as one wishes. 
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Lemma 3. Assume that sup^g^) 



Id 9hdy 



sup 



g{x,y)h{y)dy 



D 



< 1, or, more generally, that 

^ik\x-y\ 



> 0, 



9 = 9{x,y) :-- 



4:Tt\x — y\ 



Then equation (jH} has a unique solution: 

h{x) 



q{x) 



Uo[X) 



(9) 



M^) - lD9{x,y)h{y)dy' 

Formula Q holds also if condition (jH)) holds. 

Suppose that for a given h G L'^{D) condition (jH)) is not satisfied. Let us approximate 
h by an analytic function hi in D, for example, by a polynomial, so that 



!!/(/?) + 



47r 



e-'^^-''hi{x)dx\ \ < e. 



D 



Denoting hi by h again, we may assume that h is analytic in D and in a domain which 
contains D. In the following lemma we prove that it is possible to perturb h slightly so 
that for the perturbed /i, denoted /i^, condition (jS)) is satisfied, and formula Q yields a 
potential qs G L'^{D), for which inequality Q holds. 

Lemma 4. Assume that h is analytic in D and bounded in the closure of D. There exists 
a small perturbation hs of h, \\h — hs\\L2{D) < ^, such that the function 



Qs 



hs{x) 



Ug{x) - j^g{x,y)hs{y)dy 
is bounded. 

In Section 121 proofs are given. 



2 Proofs 

Proof of Lemma^ Assume the contrary. Then 

I d(3f{(3) I e-'^^-='h{x)dx = ^ WheL\D). 
Js^ Jd 

This implies 

[ f{p)e-'^^-''dl3 = Q XED. (10) 

The left-hand side of fjlO|) is the Fourier transform of a distribution /(/?)^^^^t^, where 
6{X — k) is the delta-function and A/3 is the Fourier transform variable in spherical coor- 
dinates. By the injectivity of the Fourier transform, ^TU^ implies that /(/?) = 0. □ 
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Remark 1. In this proof and in the proof of Lemma{^ below, the domain D can be taken 
such that dianiD is arbtrarily small. 



Proof of Lemma Note that Lemma |21 implies Lemma ^ It is known that 

oo 



:ii) 



i=0 



where J2e = 'Ee 



Je[r) 



f-</£+i(^), Ji{f^) is the Bessel function, Y^^^ are 



orthonormal in L {S ) spherical harmonics, Y^^mi—P) = (— 1) y^,™!/^), ct' = -, 
and the overbar stands for complex conjugate. Expand / in the Fourier series: 

oo 

f{f3) = Y.fi,mYe,miP). 

Choose L, such that 



r = \x\ 



:i2) 



Let i> L, h{x) = Y,T=o hE,mir)Yi^m{a'), a' = = \x\. 

Define hi^^ by (0). Then, by Parseval's equation, inequality Q holds. Lemma |21 is 
proved. □ 

Remark 2. Alternatively, one may look for h in the form h = 'Ej^iCjipj{x) , where 
{ipj,ipi) = 6ij and Cj = const, and minimize the left-hand side of (jU with respect to Cj, 
1 < j < J. If J is sufficiently large, the minimum will be < e. One may take not an 
orthonormal system of (fj, but just a linearly independent, complete (total) in L'^{D), 
system of functions {^j}- 

Proof of Lemma\^ The scattering solution, corresponding to a potential q, solves the 
equation 

^ika-x (13) 



u = Uq- I g{x, y)q{y)u{y)dy, Uq := e'' 

J D 

If (jHI) holds, i.e., if h corresponds to a g G L'^{D), then u = uq — Jj^ghdy. Multiply this 
equation by q and get 

q{x)u{x) = q{x)uo{x) - q{x) / g{x,y)h{y)dy. 



D 



Using (0) and solving for q, one gets provided that (jHJ holds. Condition 
||^||l2(z)) is sufficiently small. One has 



g{x,y)hdy 



D 



\x - y\ 



\h\ 



L-2(D)-, 



holds if 



(14) 
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and ^ (^/o Ji^Wj - ^'^^^^ ^ = \diamD. If, for example, -^||/i||i2(£,) < 1, then 
condition (jHJ holds, and formula © yields the corresponding potential. This explains 
the role of the "smallness" assumption. 

If /i G L'^{D) is an arbitrary function, such that condition (jSI) holds, then formula Q 
defines a potential q G L'^{D), formula © defines the function u = |, and this u solves 
equation ()13p . which is the equation for the scattering solution. Thus, the function u = ^ 
is the scattering solution, and the potential q, defined by formula corresponds to the 
given h. The above argument is valid in the case of complex-valued potential q. 

By Lemma 4, proved at the end of the paper, if h is such that condition (jHJ fails, then 
a small perturbation hs of h leads by formula Q, with hs in place of h, to a potential qs 
which generates the radiation pattern close to the desired /. 

Lemma ini is proved. □ 

Remark 3. Formula © yields, possibly, a complex-valued potential. If k > 0, then 
k"^ > is not an eigenvalue of the Schrddinger operator — + q{x) with a complex- 
valued compactly supported q G L'^{D). This follows, e.g., from the results in for 

g = o ( 4| ) as \x\ oo, and can be proved for a compactly supported q by a simple 



In the proof of Lemma|21one may take h{x) = Yl'eLo he,^mYt,m |f| j , l^^l < & < 1, where 
hi^rri do not depend on r = h{x) = for > b, and 6 > can be an arbitrary small 
number. 

Proof of Theorem^ Given an arbitrary / G L'^{S'^) and an arbitrary small e > 0, choose 
h, such that Q holds. This is possible by Lemma ^ and is done analytically in Lemma |21 
If such an h is found, then one calculates q by formula (9), provided that h is sufficiently 
small, i.e., is sufficiently small. If is arbitrary, so that conditions (jHl) or (fT^ 
are not satisfied, then q is found numerically. 

Let us explain the role of the "smallness" assumption on /. From Lemma El it follows 
that / and h are proportional, so that if 1 1/| |j;^2(-52) is sufficiently small, then 1 |/^2(£)) is 
small, and then condition (jH} is satisfied. In this case the potential is given by formula 
0. Theorem 1 is proved. □ 

Remark 4. // (jH)) fail, then formula rnay yield a q ^ LF'{D). As long as formula 
Q yields a potential q G L'^{D), p > 1, our arguments essentially remain valid. In 
our presentation we have used p = 2 because the numerical minimization in L'^-norm is 
simpler. 

The difficulty arises when formula yields a potential which is not locally integrable. 
Numerical experiments showed that this case did not occur in practice in several test 
examples in which the "smallness" condition was not satisfied. 

We prove that a suitable small perturbation hs of h in L'^{D)-norm yields by formula 
Q a bounded potential qs. This means that the "smallness" restriction on the norm of 
f can be dropped. The proof is given below. 




argument. 
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Proof of Lemma ^ Let 

N := {x : ij{x) =0,x e D}, 

where 

ij := Uo{x) - / g{x,y)hs{y)dy. 



J D 

This set is generically a hne, defined by two simultaneous equations ipj = 0,j = 1,2, 
where ijji := Keijj and 1^2 '■= Im?/'. Let 

Ns := {x : \ij{x)\ <6,x E D} 

and Ds := D \ Ns. Generically, |VV^| > c > on A^, and, therefore, by continuity, in 
A'"^. A small perturbation of h will lead to these generic assumptions. Consider the new 
coordinates 

Sl = ^1, S2 = 1p2, S3 = X3. 

Choose the origin on A^. The Jacobian f/'^^''^^''^^'', is non-singular in Ns- The vectors 



V'j/'j, j = 1,2, are linearly independent in A'^. Define hs = h in Ds and hs = in Ns- 
Let 

hs 

Qs ■■= — , 

where 

ips ■■= uo{x) - / g{x,y)hs{y)dy. 

J D 

We wish to prove that the function qs is bounded. It is sufficient to check that \iIjs\ > c > 
in A"^. By c we denote various positive constants independent of 5. One has 

\i^8\>\i^\-m>^-i{5). 

where 

I{S):=^ [ xeDs, M = max\^\. 

4?^ Jns F - y\ ^eTVi 

The proof will be completed if the estimate 

I{6) =0{S^\\n{S)\) 

is established. Let us derive this estimate. It is sufficient to check this estimate for the 
integral 

Jns \y\ Jp<5 Jo vs| + p2 

where p"^ = sf + s^. A direct calculation yields the desired estimate: / = 0(5^| ln(5)|). 
Lemma is proved. □ 
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Finally we make some remarks about ill-posedness of our algoritm for finding q given 
/. This problem is ill-posed because an arbitrary / G L'^{S'^) cannot be the scattering 
amplitude Aq{P) corresponding to a compactly supported potential q. Indeed, it is proved 
in [5], [7], that A{i3) is infinitely differentiable on 5^ and is a restriction to S"^ of a function 
analytic on the algebraic variety in C^, defined by the equation (3 ■ (3 = k"^. Finding h 
satisfying (4) is an ill-posed problem if e is small. It is similar to solving the first-kind 
Fredholm integral equation 



1 

47r 



D 



whose kernel is infinitely smooth. Our solution (7) shows the ill-posedness of the problem 
because the denominator in (7) tends to zero as £ grows. Methods for stable solutions of 
ill-posed problems (see 0) should be applied to finding h. If h is found, then q is found 
by formula (jHI), provided that ^ holds. If (jH)) does not hold, one perturbs slightly h 
according to Lemma 4, and get a potential qs by formula Q with hs in place of h. 
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